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Abstract. We examine the light-cone effect on the two- 
point correlation functions using numerical simulations 
for the first time. Specifically, we generate several sets of 
dark matter particle distributions on the light-cone up to 
z = 0.4 and z = 2 over the field-of-view of tt degree 2 from 
cosmological N-body simulations. Then we apply the se- 
lection function to the dark matter distribution according 
to the galaxy and QSO luminosity functions. Finally we 
compute the two-point correlation functions on the light- 
cone both in real and in redshift spaces using the pair- 
count estimator and compare with the theoretical pre- 
dictions. We find that the previous theoretical modeling 
for nonlinear gravitational evolution, linear and nonlinear 
redshift-distortion, and the light-cone effect including the 
selection function is in good agreement with our numerical 
results, and thus is an accurate and reliable description of 
the clustering in the universe on the light-cone. 

Key words: cosmology: theory - dark matter - large- 
scale structure of universe - galaxies: general - quasars: 
general 



1. Introduction 

In the proper understanding of on-going redshift surveys 
of galaxies and quasars, in particular the Two-degree Field 
(2dF) and the Sloan Digital Sky Survey (SDSS), it is 
essential to establish a theory of cosmological statistics 
on the light cone. This project has been undertaken in 
a series of our previous papers (Matsubara, Suto, & Sza- 
pudi 1997; Yamamoto & Suto 1999; Nishioka & Yamamoto 
1999; Suto et al. 1999; Yamamoto, Nishioka, & Suto 2000; 
Suto, Magira & Yamamoto 2000). Those papers have for- 
mulated the light-cone statistics in a rigorous manner, de- 
scribed approximations to model the clustering evolution 
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in the redshift space, and presented various predictions in 
canonical cold dark matter (CDM) universes. Their pre- 
dictions, however, have not yet been tested quantitatively, 
for instance, against numerical simulations. This is not 
surprising since it is fairly a demanding task to construct 
a reliable sample extending over the light-cone from the 
conventional simulation outputs at a specified redshift, z. 

In the present paper, we examine, for the first time, 
the validity and limitation of the above theoretical frame- 
work to describe the cosmological light-cone effect against 
the mock catalogues on the light-cone. Such catalogues 
from cosmological A-body simulations have been origi- 
nally constructed for the study of the weak lensing statis- 
tics (Hamana et al. 2000, in preparation). Applying the 
same technique (§3.1), we generate a number of differ- 
ent realizations for the light-cone samples up to z = 0.4 
and z = 2, evaluate the two-point correlation functions 
directly, and compare with the theoretical predictions. 



2. Predictions of two-point correlation functions 
on the light cone 

In order to predict quantitatively the two-point statistics 
of objects on the light cone, one must take account of 
(i) nonlinear gravitational evolution, (ii) linear redshift- 
space distortion, (iii) nonlinear redshift-space distortion, 
(iv) weighted averaging over the light-cone, (v) cosmolog- 
ical redshift-space distortion due to the geometry of the 
universe, and (vi) object-dependent clustering bias. The 
effect (v) comes from our ignorance of the correct cos- 
mological parameters, and (vi) is rather sensitive to the 
objects which one has in mind. Thus the latter two effects 
will be discussed in a separate paper, and we focus on the 
effects of (i) ~ (iv) throughout the present paper. 

Nonlinear gravitational evolution of mass density fluc- 
tuations is now well understood, at least for two-point 
statistics. In practice, we adopt an accurate fitting for- 
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mula (Peacock & Dodds 1996) for the nonlinear power 
spectrum P^[(k, z) in terms of its linear counterpart. 

Then the nonlinear power spectrum in redshift space 
is given as 



P^(k, (l) = P*{k, z)[l + 0^¥D ml [knap] 



(1) 



where k is the comoving wavenumber, and /i is the di- 
rection cosine in /c-space. The second factor in the right- 
hand-side comes from the linear redshift-space distortion 
(Kaiser 1987), and the last factor is a phenomenological 
correction for non-linear velocity effect. In the above, we 
introduce 



0(z) 



1 dlixD(z) 
b(z) dlna 



(2) 



where D(z) is the gravitational growth rate of the lin- 
ear density fluctuations, a is the cosmic scale factor, and 
the density parameter, the cosmological constant, and the 
Hubble parameter at redshift z are related to their present 
values respectively as 



Q(z) 



\(z) 



Ho 



H(z) 
H 



H(z) 



(3) 



(4) 



H{z) = iW^o(l + zf + (l - fi - A )(l + z) 2 + A (5) 

We assume that the pair-wise velocity distribution in real 
space is approximated by 



fv{vi2) 



1 



■ exp 



V2| 



"12 



(6) 



with cr P being the 1-dimcnsional pair-wise peculiar veloc- 
ity dispersion. In this case the damping term in Fourier 
space, D ve i [kjiap] , is given by 



D ve i[kna P ] = 



1 



1 + K 2 H 



-2, ,2 ' 



(7) 



where 



k(z) 



k(l + z)a r (z) 
V2H(z) 



Note that this expression is equivalent to that in Magira 
et al. (2000) but written in terms of the physical velocity 
units. 

On large scales, a r (z) can be well approximated by a 
fitting formula proposed by Mo, Jing & Borner (1997): 



n(z)H 2 (z) 

(1 + zf 



1 - 



1 + z 



D 2 {z') 



D 2 {z)J z (1 + z') 2 
dk A2 L (fc,z) 



dz' 



(9) 



We compute the pairwise velocity dispersion of par- 
ticles in iV-body simulations (see §3.1) both for SCDM 



and ACDM, whose parameters are summarized in Table 
[j], to test the accuracy of the fitting formula of the pair- 
wise velocity dispersion, eq. The measured velocity 
dispersions at z = 0, 1 and 2 are shown in Figure JTJ. The 
dotted lines in Figure [l] indicate predictions of eq. (Q) inte- 
grated over the wavenumbers existing in our iV-body sim- 
ulations. The analytical model predictions agree with our 
data within a 10% accuracy at the large separations. This 
level of agreement is as good as that found originally by 
Mo et al. (1997). Nevertheless since we are mainly inter- 
ested in the scales around l/i _1 Mpc, we adopt the follow- 
ing fitting formula throughout the analysis below which 
better approximates the small-scale dispersions in physi- 
cal units: 



, . / 740(1 + 2 )- 1 km/s 
<W) ~ j 650(1 + z)-°- 8 km/ 



for SCDM model 
km/s for ACDM model. 



(10) 



Integrating equation (Q) over [i, one obtains the 
direction-averaged power spectrum in redshift space: 



= A {k) + 2 -(3(z)B(n) + \p 2 (z)C(k) 



where 
A(k) - 
B{k) -- 
C(k) -- 



arctan(/i) 



5 
3^2 



1 - 



arctan(/i) 



3 arctan(ft) 



(11) 

(12) 
(13) 
(14) 



Adopting those approximations, the direction- 
averaged correlation functions on the light-cone are 
finally computed as 



dz ~r- \ ( i ) { z ) n o( z )\ 2 £,( x ^ z ) 

dz 



(15) 



where z m [ n and z max denote the redshift range of the sur- 



(8) vey, and 



a/ x 1 f°° nS/, ,sinfca; s , 2 „ 

£(x B ; z) = -^ P^(k, z)— k 2 dk. 
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(16) 



Throughout the present analysis, we assume a stan- 
dard Robertson - Walker metric of the form: 



ds z 



-dr 



a(t) 2 {d X 2 + S K ( X ) 2 [dd 2 



% (17) 



where Sk(x) is determined by the sign of the curvature 
K as 



(18) 



sin(V^x)/VA: [K>Q) 

Sk( X ) = { X (K = Q) 

smhiV^Kx)/^^ (K < 0) 
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Fig. 1. Pairwise peculiar velocity dispersions of dark matter particles at z = 0, 1 and 2. Dashed lines indicate the 
values predicted from the formula of MJB (eq. H), while solid lines indicate our adopted fit (eq.|10|). 



The radial comoving distance x( z ) is computed by 

x{z) = L°W) = U m (19) 

In our definition, K is not normalized to ±1 and 0, but 
rather written in terms of the scale factor at present, do, 
the Hubble constant, Hq, the density parameter, Slo and 
the dimensionless cosmological constant, Ao: 

K = a 2 H^(Q + X Q -l). (20) 

The comoving angular diameter distance D c (z) at red- 
shift z is equivalent to S (x(z)), and, in the case of 
Ao = 0, is explicitly given by Mattig's formula: 



D c (z) 



1 



a H 1 + z 1 + fl z/2 + VT+Ooi' 



(21) 



Then dV c /dz, the comoving volume element per unit solid 
angle, is explicitly given as 



dz 



S 2 K (X) 



dz 



s 2 K (x) 



ffov^l + z f + (1 - fio - A )(l + z) 2 + Ao 



(22) 



3. Evaluating two-point correlation functions 
from N-body simulation data 

3.1. Particle distribution on the light cone from N-body 
simulations 

We test the theoretical modeling against simulation re- 
sults, we focus on two spatially-flat cold dark matter mod- 



els, SCDM and ACDM, adopting a scale- invariant primor- 
dial power spectral index of n = 1 . Their cosmological pa- 
rameters are listed in Table |l|. While SCDM are known 
to have several problems in reproducing the recent ob- 
servations (e.g., de Bernardis et al. 2000), this model is 
suitable for testing the theoretical formula since the clus- 
tering evolution on the light-cone is more significant. We 
use a series of A^-body simulations originally constructed 
for the study of weak lensing statistics (Hamana et al. 
2000, in preparation). These simulations were generated 
with a vectorized PM code (Moutarde at al. 1991) modi- 
fied to run in parallel on several processors of a CRAY-98 
(Hivon 1995). They use 256 2 x 512 particles and the same 
number of force mesh in a periodic rectangular comoving 
box. We use both the small and large boxes (Table [l]) . 

The initial conditions are generated adopting the 
transfer function of Bond & Efstathiou (1984, see also 
Jenkins et al. 1998) with the shape parameter T = Qoh. 
The amplitude of the power spectrum is normalized by the 
cluster abundance (Eke, Cole & Frenk 1996; Kitayama & 
Suto 1997). 

Using the above simulation data, we generated light- 
cone samples as follows; first, we adopt a distance observer 
approximation and assume that the line-of-sight direction 
is parallel to Z-axis regardless with its (X, Y) position 
(Fig.^). Second, we periodically duplicate the simulation 
box along the Z-direction so that at a redshift z, the po- 
sition and velocity of those particles locating within an 
interval x( z ) ^ ^x( z ) are dumped, where Ax(z) is de- 
termined by the output time-interval of the original N- 
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Table 1. Parameters in iV-body simulations. 



Model 


n 


A 


h 


cr& 


Box size 


Force resolution 


Z max 












L x x L Y x £ z [fr" 3 Mpc 3 ] 


[ft -1 Mpc] 




SCDM small box 


1 





0.5 


0.6 


80 x 80 x 160 


0.31 


0.4 


SCDM large box 


i 





0.5 


0.6 


240 x 240 x 480 


0.94 


2 


ACDM small box 


0.3 


0.7 


0.7 


0.9 


120 x 120 x 240 


0.45 


0.4 


ACDM large box 


0.3 


0.7 


0.7 


0.9 


360 x 360 x 620 


1.4 


2 



Table 2. Parameter values for the polynomial evolution model of Boyle et al. (2000). 



n 


A 


Q 


P 


M B - 5 log h 


ki 


k 2 


$*[ft 3 Mpc~ 3 mag _i ] 


i 





3.45 


1.63 


-20.59 


1.31 


-0.26 


0.80 x 10 -5 


0.3 


0.7 


3.41 


1.58 


-21.14 


1.36 


-0.27 


2.88 x 10 -6 



Table 3. Summary of number of particles in each realization (real space/redshift space). 



Model Realization Total Random selection LF based selection LF based with random selection 



SCDM small box 


1 


8193106 


/ 


8216016 


10258 


/ 


10282 


125477 


/ 


125289 


8546 


/ 


8540 




2 


8291309 


/ 


8311402 


10388 


/ 


10413 


168363 


/ 


168999 


11444 


/ 


11497 




3 


8448034 


/ 


8479865 


10591 


/ 


10627 


165217 


/ 


165192 


11221 


/ 


11218 




4 


9181442 


/ 


9250736 


11533 


/ 


11618 


175769 


/ 


175773 


11927 


/ 


11927 




5 


8263119 


/ 


8324278 


10340 


/ 


10434 


178135 


/ 


177348 


12075 


/ 


12025 


SCDM large box 


1 


6253827 


/ 


6254790 


10481 


/ 


10482 


2037314 


/ 


2041146 


10348 


/ 


10363 




2 


6321816 


/ 


6319899 


10591 


/ 


10582 


2077216 


/ 


2077310 


10552 


/ 


10552 




3 


6346239 


/ 


6342617 


10626 


/ 


10623 


2090246 


/ 


2090222 


10622 


/ 


10622 




4 


6423700 


/ 


6417089 


10767 


/ 


10757 


2102122 


/ 


2099505 


10671 


/ 


10664 




5 


6298022 


/ 


6300195 


10546 


/ 


10552 


2077854 


/ 


2079776 


10553 


/ 


10564 


ACDM small box 


1 


3253963 


/ 


3224663 


7589 


/ 


7512 


43377 


/ 


42960 


8760 


/ 


8666 




2 


4326797 


/ 


4341618 


10025 


/ 


10050 


48690 


/ 


48581 


9808 


/ 


9791 




3 


4429032 


/ 


4423464 


10263 


/ 


10258 


62073 


/ 


62274 


12517 


/ 


12553 




4 


4859939 


/ 


4842481 


11245 


/ 


11201 


59105 


/ 


59022 


11903 


/ 


11890 




5 


4993640 


/ 


4988234 


11532 


/ 


11517 


72490 


/ 


72674 


14608 


/ 


14635 


ACDM large box 


1 


5358865 


/ 


5370894 


9834 


/ 


9863 


1427660 


/ 


1429062 


9588 


/ 


9593 




2 


5277031 


/ 


5286441 


9665 


/ 


9681 


1415712 


/ 


1418226 


9498 


/ 


9516 




3 


5625180 


/ 


5631157 


10322 


/ 


10326 


1507183 


/ 


1507424 


10174 


/ 


10175 




4 


5630820 


/ 


5631761 


10326 


/ 


10326 


1511963 


/ 


1511565 


10219 


/ 


10218 




5 


5606636 


/ 


5612974 


10287 


/ 


10300 


1507176 


/ 


1508490 


10174 


/ 


10187 



body simulation. Finally we extract five independent (non- 
overlapping) cone-shape samples with the angular radius 
of 1 degree (the field-of-view of ir degree 2 ), each for small 
and large boxes as illustrated in Figure ||. In this manner, 
we have generated mock data samples on the light-cone 
continuously extending up to z = 0.4 (relevant for galaxy 
samples) and z — 2.0 (relevant for QSO samples), respec- 
tively from the small and large boxes. While the above 
procedure selects the same particle at several different red- 
shifts, this does not affect our conclusion below because 
we are mainly interested in scales much below the box size 
along the Z-direction, Lz- 



In practice, we apply the above procedure separately 
in real and redshift spaces by using z Tea j and z Q b s of each 



particle (see eq. [S below) . The total numbers of particles 
in those realizations are listed in Table § 

3.2. Pair counts in real and redshift spaces 

Two-point correlation function is estimated by the con- 
ventional pair-count adopting the estimator proposed by 
Landy & Szalay (1993): 



DD{x) - 2DR(x) + RR(x) 
RR(x) 



(23) 



For this purpose, we distribute the same number of par- 
ticles over the light-cone in a completely random fashion. 
When the number of particles in a realization exceeds 10 6 , 
we randomly select 10,000 particles as center particles in 
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Fig. 2. Schematic geometry of our light-cone samples. 
The comoving distances denoted in the figure are for 
SCDM model. In the case of ACDM model, the red- 
shifts are not changed but the radial comoving dis- 
tance is 1086(3626)/i -1 Mpc and the side length is 
120(360)/i~ x Mpc for the small(large) box. 



counting the pairs. Otherwise we use all the particles in 
the pair counts. 



The comoving separation xxz of two objects located at 
z\ and zi with an angular separation 0%2 is given by 



• T i2 = x \ + x l ~ Kx\x\(l + cos 2 0i 2 ) 



Zx\X2 y 1 — Kxf y 1 — Kx\ cos Q\i 



(24) 



where x\ = D e (zi) and X2 = D c (z2). 

In redshift space, the observed redshift z Q b s for each 
object differs from the "real" one z Tea \ due to the velocity 
distortion effect: 



-^obs — ^rcal ~t~ (1 "t" ^-real)^pec : 



(25) 



where v pcc is the line of sight relative peculiar velocity 
between the object and the observer in physical units. 
Then the comoving separation s\2 of two objects in red- 
shift space is computed as 



s\ 2 = s\ + s| - Ks\s%{1 + cos 2 flia) 



-2sis 2 yi- Ksjy/1- Ks% cos 6»i2, 
where s x = D c (z ohSyl ) and s 2 = D c (z ohSy2 ). 



(26) 



3.3. Selection functions 

In properly predicting the power spectra on the light cone, 
the selection function should be specified. In this subsec- 
tion, we describe the selection functions appropriate for 
galaxies and quasars samples. 

For galaxies, we adopt a B-band luminosity function 
of the APM galaxies (Loveday et al. 1992) fitted to the 
Schechter function: 



^^*(A) Q ex P (-A) rf (A 



(27) 



1.40 x 10" 2 /i 3 Mpc 



-3 



a 



-0.97, and M* B 



with cj) 

— 19.50 + 5 log 10 h. Then the comoving number density of 
galaxies at z which are brighter than the limiting magni- 
tude Bii m is given by 



n ga i(z, < S Um ) = J 



<j){L)dL 

L(B Uxa ,z) 

= <p*T[(a + l,x(B Vim ,z)}, 



where 



x(B Vlm ,z) 



L(B Vlm ,z) 



dr(z) 



l/i- 1 Mpc 



10 



2.2-0.4Si in 



(28) 



,(29) 



and T[v, x] is the incomplete Gamma function. Figure^ 
plots the selection function defined by 



0gal(< Aim, Z) = 
with Z m in = 0.01. 



n ga i(z, < gum) 

V^mim ^ ^lim) 



(30) 




Fig. 3. Selection function of galaxies in a case of SCDM 
model on the basis of the B-band luminosity function of 
APM galaxies (Loveday et al. 1992). 



Fig. 4. Selection function of QSOs in a case of SCDM 
model on the basis of the 2dF QSO sample (Boyle et al. 
2000). 



For quasars, we adopt the B-band luminosity recently 
determined by Boyle et al. (2000) from the 2dF QSO sur- 
vey data: 



10 0.4(l-a)[M B -Af*(z)] + 10 0.4(l-/3)[M B -M*(z){ 31 ) 

In the case of the polynomial evolution model: 



M£(z) = M B (0) - 2.5(fc l2 + k 2 z 2 ), 



(32) 



and we adopt the sets of their best-fit parameters listed 
in Table | for our SCDM and ACDM. 

To compute the B-band apparent magnitude from a 
quasar of absolute magnitude Mb at z (with the luminos- 
ity distance a\{z)\ we applied the K-correction: 

B = Mb + 5 log(d L (z)/10pc) - 2.5(1 - p) log(l + z) (33) 

for the quasar energy spectrum L v oc v~ p (we use p — 0.5). 

Then the comoving number density of QSOs at z which 
are brighter than the limiting magnitude Bn m is given by 



™QSo(z, < B\ 

im ) 



M(B lh 



$(Mb, z)dM B . 



Figure [| plots the selection function defined by 

«QSo(z, < Slim) 



4>QSo{< Sii m , z) 



with z m -,„ = 0.2. 



n QSo( z min, < Aim) 



(34) 



(35) 



In practice, we adopt the galaxy selection function 
0gai(< Bn m , z) with By lm = 19 and z m ; n = 0.01 for the 
small box realizations, while the QSO selection function 
4>QSo(< Bi ilrL ,z) with Bi im = 21 and z min = 0.2 for the 
large box realizations. We do not introduce the spatial bi- 
asing between selected particles and the underlying dark 
matter, which will be discussed elsewhere. For comparison, 
we also select the similar number of particles randomly but 
independently of their rcdshifts. It should be emphasized 
here that our simulated data are constructed to match 
the shape of the above selection functions but not the am- 
plitudes of the number densities. The field-of-view of our 
simulated data, ir degree 2 , is substantially smaller than 
those of 2dF and SDSS, and we sample particles much 
more densely than the realistic number density. Since our 
main purpose of this paper is to test the reliability of 
the theoretical modeling described in section 2, and not 
to present detailed predictions, this does not change our 
conclusions below. The numbers of the selected particles 
in each realization are listed in Table § The averaged se- 
lection functions for our five realizations in real space are 
plotted as histograms in Figures || and [|. 

4. Results 

Consider first the two-point correlation functions for par- 
ticles on the light cone but without redshift-dependent se- 
lection. Figures H and^ plot those correlations for z < 0.4 
samples from small-box simulations (upper panels) and 
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Fig. 5. Mass two-point correlation functions on the light cone without redshift-dependent selection functions in SCDM 
model. Upper: z < 0.4, Lower: < z < 2.0. Left: all particles on the light cone, Right: randomly selected particles. 



for z < 2 from large-box ones (lower panels), for SCDM 
and ACDM, respectively. In these figure, we plot the av- 
erages over the five realizations (Table ||) in open circles 
(real space) and in solid triangles (redshift space), and the 
quoted error-bars represent the standard deviation among 
them. If we use all particles from simulations (left panels), 
the agreement between the theoretical predictions (solid 
lines) and simulations (symbols) is quite good. The scales 
where the simulation data in real space become smaller 
than the corresponding theoretical predictions simply re- 
flect the force resolution of the simulations listed in Table 

0- 



In order to examine the robustness of the estimates 
from the simulated data, we randomly selected N ~ 10 4 
particles from the entire light-cone volume (independently 
of their redshifts) . The resulting correlation functions are 
plotted in the right panels. It is remarkable that the es- 
timates on scales larger than ~ l/i _1 Mpc are almost the 
same. This also indicates that the error-bars in our data 
are dominated by the sample-to-sample variation among 
the different realizations. 

Next we examine the effect of selection functions. Fig- 
ures ^ and ^| plot the two-point correlation functions in 
SCDM and ACDM, respectively, taking account of the se- 
lection functions described in the subsection 3.3. It is clear 
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Fig. 6. Same as Figure | but for ACDM model. 

that the simulation results and the predictions are in good 
agreement. It should be noted that the results shown in the 
upper-left panel (intended to correspond to galaxies) have 
substantially larger error-bars compared with the corre- 
sponding ones in Figures | and |. This is an artifact to 
some extent because of the very small survey volume in 
our light-cone samples; if one applies the galaxy selection 
function which rapidly decreases as z (see, Fig.[|), the re- 
sulting structure mainly probes the universe at z < 0.1 
and thus large-scale nonlinearity or variation for the differ- 
ent line-of-sight becomes significant. If we are able to use 
the same number of particles but extending over the much 
larger volume, the sample-to-sample variations should be 
substantially smaller. This interpretation is supported by 



x [h _1 Mpc] (comoving) 



the upper-right panel where we randomly sample N ~ 10 4 
particles from those used in the upper-left panel. Despite 
the fact that the number of particles is only 5% (20%) of 
the original one for SCDM (ACDM) model, the resulting 
correlation functions and their error-bars remain almost 
unchanged. The lower panels corresponding to QSOs show 
the similar trend. 



5. Conclusions and discussion 

We have presented detailed comparison between the the- 
oretical modeling and the direct numerical results of the 
two-point correlation functions on the light-cone. In short, 
we have quantitatively shown that the previous theo- 
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Fig. 7. Mass two-point correlation functions on the light cone for particles with redshift-dependent selection functions 
in SCDM model. Upper: z < 0.4, Lower: 0.2 < z < 2.0. Left: with selection function whose shape is the same as that 
of the B-band magnitude limit of 19 for galaxies (upper) and 21 for QSOs (lower). Right: randomly selected N ~ 10 4 
particles from the particles in the left results. 



retical models by Yamamoto & Suto (1999) and Ya- 
mamoto, Nishioka & Suto (1999) are quite accurate on 
scales l/i _1 Mpc < x < 20/i~ 1 Mpc where the numerical 
simulations are reliable. It is also encouraging that this 
conclusion remains true even for the particle number of 
around 10 4 . In fact, the error-bars in our estimates of 
the two-point correlation functions are dominated by the 
sample-to-sample variance due to the limited angular-size 
(7r-degree 2 ) and thus the limited volume. 

In order for the more realistic evaluation of the statis- 
tical and systematic uncertainties, one needs mock light- 



cone data samples with a much wider sky coverage. More 
importantly such datasets enable one to access the effect of 
biasing on the two-point correlation functions on the light- 
cone. Since our present study indicated that all the physi- 
cal effects except for the biasing are well described by the 
existing theoretical models, it is very interesting to exam- 
ine in detail how to extract the effect of the galaxy/QSO 
biasing from the upcoming redshift survey on the basis of 
the above mock samples. We plan to come back to these 
issues with larger simulation datasets in near future. 
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